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Abstract 

Neutralinos annihilating in the center of the Sun or the Earth may give rise to a 
detectable signal of neutrinos. We derive the indirect detection rates for neutrino 
telescopes in the minimal super symmetric extension of the standard model. We 
show that even after imposing all phenomenological and experimental constraints 
that make the theories viable, regions of parameter space exist which can already 
be probed by existing neutrino telescopes. We compare with the discovery po- 
tential of supersymmetry at LEP2 as well as direct detections and point out the 
complementarity of the methods. 

1 Introduction 

Supersymmetric neutralinos with masses in the GeV-TeV range are among the lead- 
ing non-baryonic candidates for the dark matter in our galactic halo. One of the most 
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promising methods for the discovery of neutrahnos in the halo is via observation of en- 
ergetic neutrinos from their annihilation in the Sun and/or the Earth 0. Through 
elastic scattering with the atomic nuclei in the Sun or the Earth, a neutralino from the 
halo can lose enough energy to remain gravitationally trapped 0. Trapped neutralinos 
sink to the core of the Sun or the Earth where they annihilate into ordinary particles: 
leptons, quarks, gluons and - depending on the masses - Higgs and gauge bosons. Be- 
cause of absorption in the solar or terrestrial medium, only neutrinos are capable of 
escaping to the surface. Neutralinos do not annihilate into neutrinos directly [Q, but 
energetic neutrinos may be produced via hadronization and/or decay of the direct anni- 
hilation products. These energetic neutrinos may be discovered by terrestrial neutrino 
detectors. 

In this letter, we consider Cerenkov neutrino telescopes. They consist of large un- 
derground arrays of photo-multipliers to detect the Cerenkov light emitted by muons 
generated in charged- current interactions of neutrinos with the medium surrounding the 
detector. Underground Cerenkov detectors, originally built to search for proton decay, 
have already started to explore (and constrain) neutralinos as dark matter candidates 

p. A new generation of much larger neutrino telescopes, utilizing large volumes of 
natural water or ice (for a review, see 0) is currently under construction, which will 
increase the sensitivity for high-energy neutrino point sources by an order of magnitude 
or more. 

The prediction of muon rates is quite involved: we compute neutralino capture 
rates in the Sun and the Earth, fragmentation functions in basic annihilation processes, 
propagation through the solar or terrestrial medium, charged current cross sections and 
muon propagation in the rock, ice or water surrounding the detector. In addition, there 
may be scattering of the Cerenkov photons generated by the muons, due to impurities 
in the medium. 

2 Definition of the Super symmetric Model 

The minimal = 1 supersymmetric extension of the standard model is defined by the 
following superpotential and soft supersymmetry-breaking potential (for notation and 
details, see Ref. [||) 

W = 6., (-e*j,YE\Hi - d*j,Ync^LHl + UrYuc^lM - f^Him) , (1) 

-BiiH{Hi + h.c.) 

+H{*m\H{ + W*mlHi 
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M2 ( W^W^ + 2W+W' 



(2) 



Here i and j are SU(2) indices (ei2 = +1), Y's, A's and M's are 3x3 matrices in 
generation space, and the other boldface letters are vectors in generation space. 

Electroweak symmetry breaking is caused by both H\ and if| acquiring vacuum 
expectation values, 



(^2) = ^2, 



(3) 



with g'^iv'l + V2) = 2m^, with the further assumption that vacuum expectation values 
of all other scalar fields (in particular, squark and sleptons) vanish. This avoids color 
and/or charge breaking vacua. 

To reduce the number of free parameters, we make simplifying unification assump- 
tions for the gaugino mass parameters 



Ml 
M2 



tan^ 6'^M2, 



sin^ 9,nOt 



(4) 
(5) 



and a simple ansatz for the soft supersymmetry-breaking parameters in the sfermion 
sector: 



Ad 
As 



diag(0,0,At), 
diag(0,0, A), 
0. 



M, 



M, 



mol. 



(6) 
(7) 
(8) 
(9) 



We choose as independent parameters the mass ttla of the CP-odd Higgs boson, the 
ratio of Higgs vevs tan/5 = V2/V1, the gaugino mass parameter M2, the higgs(ino) mass 
parameter /i, and the quantities mo. At and A}, above. 

We include one-loop radiative corrections and two-loop leading-log contributions to 
the Higgs mass matrices using the effective potential approach described in 0. 

The neutralinos Xi* are linear combinations of the neutral gauginos 5, and of the 
neutral higgsinos ifg- ^^is basis, their mass matrix 



Al,2,3,4 
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(10) 



can be diagonalized analytically to give four neutral Majorana states. 



(11) 
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the lightest of which, to be called is then the candidate for the particle making up 
(at least some of) the dark matter in the universe. 

We consider only models that satisfy all accelerator constraints on supersymmetric 
particles and couplings; in particular, the measurement of the b — ^ S7 process at the 



Cornell accelerator |I0[, which provides important bounds. However, to show the impact 
of present dark matter searches, we plot also models which are excluded by these non- 
accelerator searches. 



3 Relic Density and Capture Rates 

For each model allowed by the accelerator constraints we calculate the relic density of 



neutralinos Vt^h . We use the formalism in Ref. |11] to carefully treat resonant anni 



hilations and threshold effects, keeping finite widths of unstable particles, including all 
two-body annihilation channels of neutralinos. The annihilation cross sections were de- 
rived using a novel helicity projection technique fl^, and were checked against published 
results for several of the subprocesses. 

In this letter, we keep only models in which the neutralino density does not overdose 
the Universe and in which neutralinos can make up the totality of the galactic dark 
matter. Namely we require f^gai.DM^^ < f^^^^ ^ where (somewhat arbitrarily) we 
choose figai.DM^^ = 0.025. We have adopted a local dark matter density of 0.3 GeV/cm^. 

The capture rate in the Earth is dominated by scalar interactions, and presents 
kinematic enhancements whenever the mass of the neutralino almost matches one of 
the heavy elements in the Earth. For the Sun, both axial interactions with hydrogen 
and scalar interactions with heavier elements are important. For both the Sun and the 
Earth we use the convenient approximations available in |13|. 



4 Muon fluxes from neutralino annihilations 

Neutralinos in the core of the Sun and/or Earth can annihilate to a fermion-antifermion 
pair, to gauge bosons, Higgs bosons and gluons (xx ~^ qq, gg, qqg, W~^W~ , 

Z^Z^, Z^H^, W^H^, H^H^). These annihilation products will hadronize and/or decay, 
eventually producing high energy muon neutrinos. Since the rate of muons in a neutrino 
telescope is approximately proportional to the neutrino energy squared (since both the 
cross section and the muon range are approximately proportional to the energy), the 
annihilation channels with the hardest neutrino spectra will be the most important, i.e. 
W~^W~ , Z'^Z^, ti, etc. In our calculation of the neutrino fluxes we have however included 
all annihilation channels (except gluons since they give very soft neutrino spectra). 

With Monte Carlo simulations we have considered the whole chain of processes from 
the annihilation products in the core of the Sun or the Earth to detectable muons at 
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the surface of the Earth [|n 

the hadronization and decay of the annihilation products using Jetset 7.4 
neutrino interactions on their way out of the Sun and of the charged-current neutrino 
interactions near the detector using Pythia 5.7 [1^, and finally of the multiple Coulomb 



We have performed a full Monte Carlo simulation of 

H, of the 



scattering of the muon on its way to the detector using distributions from Ref. 

PI 
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With respect to calculations using Ref. (e.g. Ref. [|T8|),this Monte Carlo treat- 
ment of the neutrino propagation through the Sun bypasses simplifying assumptions 
previously made, namely neutral currents are no more assumed to be much weaker than 
charged currents and energy loss is no more considered continuous. For details on this 
treatment, see Ref. [113, 



m 



The muon flux at a detector has been simulated for a set of neutralino masses (m^ = 
10, 25, 50, 80.3, 91.2, 100, 150, 176, 200, 250^ 350, 500, 750, 1000, 1500, 2000, 3000 and 
5000 GeV) and annihilation channels (cc, bb, tt, t'^t" , W~^W~ and Z^Z^). For each 
mass and channel, 2.5 x 10^ annihilations have been simulated. For masses other than 
those simulated, an interpolation is performed and the muon flux from channels other 
than those listed above are easily calculated since all other annihilation products decay 
to these particles (lighter quarks, electrons and muons do not contribute significantly to 
the neutrino flux). For the Higgs bosons, which decay in flight, an integration over the 
angle of the decay products with respect to the direction of the Higgs boson is performed. 
Given the branching ratios for different annihilation channels it is then straightforward 
to obtain the muon flux above any given energy threshold and within any angular region 
around the Sun or the center of the Earth. 



5 Indirect detection rates 



To illustrate the potential of neutrino telescopes for discovery of dark matter through 
neutrinos from the Earth or the Sun, we present the results of our full calculation. We 
show together results obtained with one 'normal' scan in the parameter space letting 
/i, M2, tan/5, rriA, mo, Ab and At vary at random between generous bounds and one 
'special' scan where we have been more restrictive on the A mass. 



/i G [-5000, 5000] GeV 
M2 G [-5000, 5000] GeV 
tan/? G [1.2,50] 
tha G [0, 1000] GeV 
mo G [100, 3000] GeV 
Ab G [-3,3]mo 
At G [-3, 3]mo 



'normal' 



/i G [-5000, 5000] GeV 
M2 G [-5000, 5000] GeV 
tan/? G [1.2,50] 
ruA G [0, 150] GeV 
mo G [100, 3000] GeV 
Ab G [-3,3]mo 
At G [-3,3]mo 



'special' 



We recall that the density of points in the figures reflects our choices for scanning the 
parameter space, and is therefore subjective (for a discussion on this see Ref. M). 
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In Fig. |1| we show our predictions for the indirect detection rates as a function of 
neutrahno mass. The horizontal hues are the best present hmits for indirect searches 
and come from the Baksan detector [Q. The hmits are (^^"-"'^^ < 2.1 x 10~^^ cm~^ 
s-i and $2"" < 3.5 x IQ-^^ cm ^ s ^ at 90% confidence level and integrated over a 
half-angle aperture of 30° with a muon energy threshold of 1 GeV. Quite a few models 
with high neutrino rates are already ruled out by the Baksan data. These models have 
large tan (3 > 30, low H2 mass, and large mixing for stop squarks. In Fig. |l| it can 
also be seen that a neutrino telescope of an area around 1 km^, which is a size currently 
being discussed, would have a large discovery potential for supersymmetric dark matter. 
We remind that, after subtraction of the atmospheric neutrino fluxes (by means of an 
on-source off-source difference, e.g.), the remaining background is due to high-energy 
neutrinos produced by cosmic ray collisions in the solar atmosphere, and is at the level 
of 15 muons (> 1 GeV)/kmVyr pg. 

In Fig. we show the muon rates versus Q^h"^. The general trend is that large ^x^"^ 
corresponds to lower rates, as would be expected from crossing symmetry between anni- 
hilation and scattering cross sections. Note, however, the large spread of the predicted 
rates for a given value of Q^h"^. 

In Fig. 1^ a comparison is made between the predicted rates from the Earth and from 
the Sun for the same set of models. In the region yet to explore, the rate from the Sun 
is generally larger than the rate from the Earth. 

In Fig. ^the indirect detection rate is compared to the direct detection rate in ^^Ge. 
As can be seen, there is a correlation between the two, although for the Sun it is not as 
strong as for the Earth, where a high capture rate is due to a large scalar cross section, 
which also means a high rate in Germanium^. Without forgetting the huge spread, we 
see that for a given factor of improvement in sensitivity, indirect detection from the Sun 
generally scores better than direct detection, which in turn generally scores better than 
indirect detection from the Earth. 

Note that the muon rates in real experiments may be significantly lower (by as much 
as an order of magnitude for neutralinos in the lower mass range) due to the need to 
impose a higher energy threshold for the signal than assumed here. We have taken 1 
GeV for neutrino telescopes which is true for a small scale detector like Baksan; for a 
kilometer-scale array it is more likely to be tens of GeV. Likewise, the Germanium rate 
given is the integrated rate from zero recoil energy to the kinematical limit. Present-day 
detectors typically only sample a small range of recoil energies. 

Since the special scan employs fairly low values of the A mass, one may wonder 
whether the A or the lightest Higgs boson H2 would be light enough to be discovered 
or excluded at LEP2. In particular, one could expect this for the models that give 

^Recently, a paper appeared (V.A. Bednyakov et al., bep-ph/960626l[ ), where no high Germanium 
rates were found in a variant of the special scan employed here. This seems to be due to their use of a 
more restricted model, imposing unification conditions of the scalar mass parameters at the GUT scale, 
which we do not. 



6 



high neutrino rates, since a low Higgs mass gives a large spin-independent scattering 
cross section of (mixed) neutralinos and therefore large capture and annihilation rates. 
Using the expected exclusion limits for LEP2 given in Ref. ||2T]] we find, indeed, that the 
models in our sample with the highest rates are within reach of LEP2, because of the 
large cross section for e+e~ — > AH2. Technically this is due to the fact that these models 
have a large value of the factor cos^(/? — a) that governs the ZAH2 coupling (with a 
being the CP-even Higgs mixing angle). In Fig. ^ we show the indirect detection rates 
versus 771^2- the models that can be probed by LEP2 are shown in the upper part of 
the figure, the others in the lower part. We have used the combination of all four LEP 
experiments and assumed a total integrated luminosity of 150 pb~^ at 192 GeV. We 
have also included the possible LEP2 lower limit of 95 GeV on the mass of the chargino 
and the bounds coming from the e~^e^ ZH2 Zbb process. LEP2 will probe all 
models in our scan with < 90 GeV, and in particular most of the models already 
ruled out by Baksan. On the other hand, there are many models giving quite large 
indirect detection rates that cannot be probed by LEP2, but that would be accessible 
at a large neutrino telescope. 

In conclusion, indirect dark matter searches and LEP2 probe complementary re- 
gions of the supersymmetric parameter space. Moreover, direct detection (see |§]) is 
reaching a sensitivity that allows some models to be excluded, with somewhat different 
characteristics than those probed by the other methods. This illustrates a nice com- 
plementarity between direct detection, indirect detection and accelerator methods to 
bound or confirm the minimal supersymmetric standard model. 
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Figure 1: The indirect detection rates from neutralino annihilations in a) the Earth and 
b) the Sun versus the neutralino mass. The horizontal line is the Baksan limit 0. Filled 
points are from a 'normal' scan and open points are from a 'special' scan with < 150 
GeV. 
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Figure 2: The indirect detection rates from neutralino annihilations in a) the Earth and 
b) the Sun versus Q^h^. The horizontal line is the Baksan limit 0. Filled points are 
from a 'normal' scan and open points are from a 'special' scan with < 150 GeV. 
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Figure 3: The indirect detection rates in the Earth versus the indirect detection rates 
in the Sun. The horizontal and vertical lines are the Baksan limits The dashed line, 
indicating equal rates, is shown just for convenience. Filled points are from a 'normal' 
scan and open points are from a 'special' scan with < 150 GeV. 
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Figure 4: The indirect detection rates from neutralino annihilations in a) the Earth and 
b) the Sun versus the direct detection rate in ^^Ge. The horizontal line is the Baksan 
limit P|. Filled points are from a 'normal' scan and open points are from a 'special' 
scan with < 150 GeV. 
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Figure 5: The indirect detection rates from neutralino annihilations in a) the Earth and 
b) the Sun versus mj^j. The horizontal line is the Baksan limit 0. The upper part 
are models which can be excluded at LEP2 after 150 pb~^ of running at 192 GeV and 
the lower part are models which cannot be excluded by LEP2. The open triangles are 
models that can be excluded due to no Higgs discovery and the open squares are models 
that can be excluded due to no chargino discovery. 
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